Spin- wave spectrum in La2Cu04 
double occupancy and competing interaction effects 
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The recently observed spin- wave energy dispersion along the AF zone boundary in La2CuC>4 
is discussed in terms of double occupancy and competing interaction effects in the t — t' Hubbard 
model on a square lattice. 
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I. INTRODUCTION 

Recently high resolution inelastic neutron scattering 
studies of the spin-wave spectrum have been carried out 
in the square-lattice spin-1/2 antiferromagnet La2Cu04. 
Spin-wave energies were obtained along high symme- 
try directions in the two-dimensional (2D) Brillouin zone. 
With wavevectors identified by their coordinates (h, k) in 
the 2D reciprocal space of the square lattice, these stud- 
ies reveal a noticeable spin-wave dispersion along the AF 
zone boundary, with a maximum of 335 meV near the 
point X having Q = (1/2, 0) and a minimum of 285 meV 
near Q = (3/4, 1/4) at T = 10 K. 

The spin-wave spectrum for the Quantum Heisenberg 
antiferromagnet (QHAF) with only nearest-neighbour 
(NN) exchange coupling J, conventionally obtained us- 
ing the linear spin- wave theory, shows no dispersion along 
the AF zone boundary. The simplest explanation for the 
observed spin-wave dispersion involves a next-nearest- 
neighbour (NNN) ferromagnetic exchange coupling J' be- 
tween Cu spins, and best fit was found with J = 104. 1±4 
meV and J' = -18 ± 3 meV at T = 10 K. § 

It was pointed out Q that a more natural explanation 
of the data can be obtained in terms of the one-band 
Hubbard model, for which the strong coupling expansion 
in powers of the NN hopping term t up to fourth order 
0(t 4 ) yields an effective spin Hamiltonian with antiferro- 
magnetic exchange interactions between NN, NNN, and 
NNNN pairs of spins, as well as a cyclic ring exchange 
interaction coupling four spins at the corners of a square 
plaquette. |^| |J The quantitative description within lin- 
ear spin- wave theory of both the spin- wave energies and 
intensities was taken as a strong demonstration that the 
one-band Hubbard model is an excellent starting point 
for describing the magnetic interactions in cuprates. Fits 
were found to be indistinguishable from those for vari- 
ables J and J', and yielded t = 0.30 ± 0.02 eV and 
U = 2.2 ± 0.4 eV at T = 10 K, so that the characteristic 
ratio U/t — 7.3. Recently, excellent fit has also been ob- 
tained with the spin-wave energy for the square-lattice 
Hubbard model, evaluated numerically in the Random 
Phase Approximation (RPA). || 

The one-band Hubbard model with only NN hopping 
possesses particle-hole symmetry, implying identical be- 
haviour for hole and electron doping. However, there 



exist significant differences in the magnetic properties 
of hole-doped and electron-doped cuprates. In partic- 
ular, AF order persists in the electron-doped cuprate 
Nd2- x Ce x Cu04 up to a doping concentration of about 
15%, HQ whereas only 2% hole concentration destroys 
AF order in La2- x Sr x Cu04. Thus, the Hubbard model 
with only NN hopping is unable to describe the magnetic 
properties of doped cuprates. 

Recently, the magnetic phase diagram of the t—t' Hub- 
bard model has been obtained in the t' — U space, || for 
electron and hole doping and both signs of the NNN hop- 
ping term t' . The phase diagram shows that for an ap- 
propriate sign of t' , the AF state survives electron doping 
(up to x « 20% for |i'/t| = 0.25 and U/t w 8), whereas 
for any finite hole doping the AF state is destroyed. Thus 
the t — t' Hubbard model appears to be the simplest cor- 
related electron model capable of describing the magnetic 
correlations in cuprates, even at the doped level. 

The need for more realistic microscopic models, which 
include NNN hopping etc., has also been acknowledged 
recently from band structure studies, photoemission data 
and neutron-scattering measurements of high-T c and re- 
lated materials. [p|~[T2"f Estimates for \t' /t\ range from 0.15 
to 0.5. 

Now, the NNN hopping term t' generates antiferromag- 
netic NNN interaction of comparable magnitude (t 2 /t 2 is 
of the order of J'/ J). The resulting competition between 
the exchange interactions will modify the values of t and 
U obtained earlier. [Q In this brief report we obtain the 
modified one-band Hubbard model parameters by fitting 
the spin- wave spectrum for La2Cu04 with that for the 
t — t' Hubbard model. We have used the RPA approach 
to directly obtain the spin-wave propagator in the AF 
state of the t — t' Hubbard model. This approach has 
the advantage of being applicable in the intermediate- 
and weak-coupling regimes as well, as it eliminates the 
need for going through the strong coupling expansion to 
obtain the effective spin model. 

Finite-?/ effects of double occupancy on the spin-wave 
dispersion in the half-filled Hubbard model on a square 
lattice have been studied earlier at the RPA level, both 
analytically by systematically expanding in powers of 
t/U, |13[ a nd also numerically through an exact compu- 
tation - WM This latter approach, desirable in view of the 
slow, asymptotic convergence of the t/U expansion, was 
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employed to study the spin- wave properties in the weak 
coupling limit, where it was found that the spin- wave en- 
ergy scale gets compressed by the decreasing charge gap. 

In ref. [13], the spin- wave propagator at the RPA level, 
X~ + (q,^) = X°(q,^)/l - Ux°{q.,w), was obtained by 
keeping all terms up to 0(t 4 /U 4 ) (and some terms even 
up to 0(t 6 /U 6 )) in the zeroth-order particle-hole propa- 
gator X (Qj uj )- The terms causing spin- wave dispersion 
along the AF zone boundary are of order t 2 /U 2 , showing 
the importance of charge fluctuations, as the double oc- 
cupancy term (n^ n^) is of the same order. This RPA 
approach is isomorphous to the linear spin-wave analysis 
P,fl5| of the effective spin- 1/2 Heisenberg model obtained 
from the strong-coupling expansion up to order t 4 , where 
the higher order exchange terms arise from the coherent 
motion of electrons beyond the NN sites. 



II. SPIN- WAVE SPECTRUM 

Including the contribution in x°(q, u>) of the NNN hop- 
ping term t' up to the leading order (t 2 /U 2 ) only, fl^ ] 
and neglecting terms of order t 2 t 2 /U 4 , and t 4 /U 4 which 
are an order of magnitude smaller than t 4 /U 4 , we obtain 
the spin-wave propagator 
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and the spin-wave energy 
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where 7 q = (cos q x + cos q y ) /2 and 7 q = cos q x cosq y , and 
M = 1 — 8t 2 /f/ 2 is the HF-level sublattice magnetization. 
Here q = Q— (n, w) is the wavevector measured from the 
zone center. The exchange energy scale J = 4t 2 /U, and 
the matrix elements a and b in Eq. (1) are given by 
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The form of Eq. (1) ensures that the spin susceptibility 
sum rule, / ^[x + "(w)-x" + (w)] = Ma z , is obeyed. g| 
This feature of the spin-wave propagator was also noted 
when quantum corrections were included, Jl8| where it is 
the one-loop-level sublattice magnetization, reduced by 
quantum spin fluctuations, which appears as the overall 
factor. 




FIG. 1. The \t'/t\ value below which w(tt,0) > w(n/2,n/2), 
as seen for La2Cu04 in the neutron-scattering studies. 

In Eq. (2) for the spin-wave energy, the double- 
occupancy term of order t 2 /U 2 and the competing- 
interaction term of order t 2 /t 2 are of similar magnitude. 
While both terms remove the degeneracy in the spin- 
wave spectrum along the AF zone boundary (7 q = 0), 
they give rise to opposite dispersion along the AF zone 
boundary. Thus, while double occupancy leads to max- 
imum and minimum spin-wave energies at points q = 
(±7r, 0), (0, ±7r) and (±vr/2, ±tt/2) respectively, the com- 
peting interaction favors an exactly opposite trend. 



III. DETERMINATION OF t and U 

From Eq. (2), the spin- wave velocity c = 
limq^olujq/q), and the spin- wave energies at points (ir, 0) 
and (7r/2,7r/2) are obtained as: 
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A momentum-independent multiplicative renormaliza- 
tion factor Z c — 1.18, arising from quantum spin fluc- 
tuations, (l^] has been included. Equations (4-6) clearly 
show that for U/t ~ 8 and t'/t ~ 1/4 (as appropriate 
for the cuprates), the double occupancy and compet- 
ing interaction terms are of the same order of magni- 
tude. It is also interesting to note that c and lj(tt/2, tt/2) 
are directly proportional to each other. Indeed, the 
spin- wave spectrum for La2CuC>4 shows that both c and 
a; (7r/2,7r/2) are nearly unchanged with temperature. |lj] 
Interestingly, from Eqs. (4) and (6) we find that c and 
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lu(tt/2,tt/2) have a vanishing derivative with respect to 
t/U (which should change with temperature) at U/t ~ 6, 
which is very close to the U/t value obtained from our 
fit. 

As the maximum in the spin- wave energy spectrum of 
La2CuC>4 is clearly seen to occur at (n,0), the double 
occupancy term is actually dominant, showing the im- 
portance of charge fluctuations in cuprates. From this 
information, an upper limit for the ratio \t'/t\ can be de- 
duced by comparing the spin-wave energies uj(it,0) and 
ui(ir/2, it/2). From Eqs. (5) and (6), we obtain the con- 
dition t 2 /t 2 < 6t 2 /U 2 , which is shown in Fig. 1. For 
U/t ~ 8, the ratio \t'/t\ must be smaller than ~ 0.3. 

Taking the key values of the maximum spin-wave en- 
ergy 0) and the spin-wave velocity c from the spin- 
wave spectrum, we solve for t/U and J from Eqs. (4) and 
(5), and thus obtain t and U. With u)(tt, 0) = 335 meV 
and c = 200 meV at T = 10 K, and taking \t'/t\ = 0.25, 
we obtain t 2 /U 2 = .022, U/t = 6.7, t = 0.34 eV, and 
U = 2.3 eV. Figure 2 shows the spin-wave spectrum cal- 
culated from Eq. (2) with these t and U values, and the 
experimental points for La2CuC>4 at 10 K. Without t', 
the corresponding values are t 2 /U 2 = .019, U/t = 7.3, 
t = 0.29 eV, and U — 2.2 eV in agreement with the re- 
sults of ref . [1] , and showing a nearly 20% enhancement 
in the double occupancy factor t 2 /U 2 due to the com- 
peting effect of if. With u(ir, 0) = 320 meV and c = 200 
meV at T = 295 K, we obtain U/t = 7.1, t = 0.34 eV, 
and U = 2.4 eV. 

In the approach adopted in ref. [1], the effective spin 
5 = 1/2 quantum Heiscnberg Hamiltonian obtained by 
the strong-coupling expansion of the Hubbard model con- 
tains all physical processes up to the order t 4 level, and 
the subsequent use of the linear spin- wave approximation 
to obtain the spin-wave spectrum neglects the quantum 
corrections of order 1/S. In our approach these two steps 
are carried out in exactly reverse order. The evaluation 
of the spin- wave propagator in the random phase approx- 
imation neglects the quantum corrections (of equivalent 
order 1/fif = 1/2S within the inverse-degeneracy expan- 
sion of the generalized TV-orbital Hubbard model) [Q . 
The evaluation of X (q, w ) is next carried out, either an- 
alytically via a systematic expansion in powers of t/U, 
or numerically through an exact computation. 

The spin-wave dispersion along the AF zone boundary, 
caused by the double occupancy effect at finite U and 
partially neutralized by the competing interaction effect 
of t' , results in a broadening of the spin- wave spectrum 
at the high-energy end, resembling a spin- wave damp- 
ing effect. Implication of this broadening to the two- 
magnon Raman scattering in La2Cu04 [^,^l| has been 
discussed earlier. Jl^JI^] The importance of charge fluctu- 
ation effects (in terms of the cyclic exchange terms) have 
also been pointed out in infrared absorption experiments, 
p3| , p4j and magnetic properties of the related compound 

Sl'i4Cu2404i. |2a] 
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FIG. 2. The calculated spin-wave energy spectrum o) q from 
Eq. (2) (line) and spin-wave energies for La2Cu04 at T = 10 
K (triangles) from ref. [1]. 

In conclusion, the competing interaction effect of the 
NNN hopping term t', which renders the Hubbard model 
more realistic for describing the magnetic properties of 
the doped cuprates, partially neutralizes the spin- wave 
dispersion along the AF zone boundary caused by the 
double occupancy effect. Due to this competition the 
double occupancy factor t 2 /U 2 must be significantly en- 
hanced by nearly 20% in order to fit the observed spin- 
wave spectrum in La2Cu04. 
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